Controlled coupling of counterpropagating whispering-gallery modes by a single 
Rayleigh scatterer: a classical problem in a quantum optical light 
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We present experiments where a single subwavelength scatterer is used to examine and control 
the back-scattering induced coupling between counterpropagating high-Q modes of a microsphere 
resonator. Our measurements reveal the standing wave character of the resulting symmetric and 
antisymmetric eigenmodes, their unbalanced intensity distributions, and the coherent nature of their 
coupling. We discuss our findings and the underlying classical physics in the framework common to 
quantum optics and provide a particularly intuitive explanation of the central processes. 



The radiative properties of atoms can be strongly mod- 
ified by coupling them to resonators [ij. A historical 
corner stone of this field of research, known as Gavity 
Quantum Electrodynamics (CQED), was set in 1946 by 
E. M. Pureed who proposed that the radiation rate of 
an oscillating dipole at wavelength A can be enhanced 
by a factor F — 'iQX'^ / An^Vm in a resonant cavity of 
quality factor Q and mode volume Vm 13 • This so- 
called Pureed effect holds in the dissipative weak cou- 
pling regime where the cavity finesse is small so that the 
atomic radiation remains dominated by its coupling to 
the bath of the electromagnetic modes. In the strong 
coupling regime, coherent exchange of energy between 
the atom and the resonator causes the atomic resonance 
to lose its identity and to become replaced by a doublet. 
These phenomena have been studied for more than three 
decades 0, H, 0, [E| although the in-situ manipulation of 
a single emitter in a single mode of a high-Q microres- 
onator remains a challenge 0, 0|. In this Letter, we con- 
sider the controlled coupling of a classical nano-object to 
a high-finesse whispering-gallery mode (WGM) microres- 
onator. We discuss both theoretically and experimentally 
the resulting coherent coupling between two degenerate 
counterpropagating WGMs and the modification of the 
Rayleigh scattering rate. Our findings show that the con- 
cepts of the strong and weak coupling play a central role 
even in this fully classical system. 
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FIG. 1: a) Whispering gallery modes of a microsphere are 
excited via a prism. A glass fiber tip can be positioned in 
(r, 6, (f)) close to the sphere surface. The resonator spectrum 
can be recorded on PDl in transmission, on PD2 from global 
scattering out of the sphere via a multimode fiber, and on 
PD3 through the fiber tip. b) An example of the doublet 
spectrum recorded on PD2. 

The resonators in our work consist of microspheres 
melted at the end of silica fibers Q . Such spheres support 
very tiigh-Q WGMs and have been studied by several 
groups |3, [1, Q. About ten years ago, it was discovered 
that the high-Q resonances of these cavities are often 
composed of doublets Such a mode splitting has 

been since discussed in conjunction with various WGM 
resonators [1, [111, [H, It turns out that mode 

splitting has been observed in other ring resonators and 
has been explained as the result of the coupling between 
the electric fields Ec and Ecc of the degenerate clockwise 
(c) and counter clockwise (cc) modes via back scattering. 
The new superpositions states (-I-) and (— ) are described 
by 



E+ = aE, + bE,, ■ E^ = aEc - bE,. 



(1) 
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Here a and b are complex coefficients. In the simplest 
case, the coupling between Ec and Ecc can be caused 
by a refiector [ij, In the case of WGM resonators, 
however, it has been suggested that backscattering from 
a distribution of residual subwavelength inhomogeneities 
in the glass matrix or on its surface is the source of 
this coupHng (1, i, [TD]. The orders of magnitude of 
the doublet splitting can be correctly estimated from 
classical electrodynamic considerations following this hy- 
pothesis [i, i, [Tfl, [H, E m. Nevertheless, the direct 
link between the spectral features of a doublet and the 
nanoscopic details of the backscattering sources has not 
been demonstrated experimentally, and a proper treat- 
ment of the losses inflicted by the scatterers is missing 
in the literature. In fact, an intuitively perplexing and 
interesting question arises in this context: given that the 
radiation of a subwavelength scatterer is nearly isotropic 
and that the angle subtended by a typical cavity mode is 
merely about 10~^ rad ^17l|, how could the rate of scat- 
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tering back into a cavity mode dominate the rate of scat- 
tering out of the resonator to ensure the population of 

E+ and E^l 

The schematics of our experimental arrangement is 
shown in Fig. [Iji. A narrow-band diode laser (A = 
670 nm, linewidth < 300 kHz, tuning range « 60 GHz) 
was used to excite the WGMs via a prism. Photode- 
tector PDl was used to record spectra in transmission 
through the prism whereas PD2 captured the light glob- 
ally scattered out of the microsphere into a multimode 
fiber. Figure[T)3 displays a typical doublet with a splitting 
of 29 MHz and Q ~ 8 X 10^ recorded on PD2. The peaks 
represent the intensities and l-E-P of the symmet- 

ric and antisymmetric eigenmodes. Following the proce- 
dure described in Ref. [l8| , we applied scanning near-field 
optical microscopy (SNOM) to map the spatial intensity 
distribution of the WGMs on PD3 and to identify the 
fundamental mode of the resonator, which exhibits a sin- 
gle intensity maximum in the r and 6 directions. In our 
previous works, we have shown that an uncoated glass 
tip might broaden and shift cavity modes depending on 
their Q and on the tip size 19]. Here we demonstrate that 
a subwavelength tip can modify or induce the coupling 
between the degenerate c and cc microsphere modes. 

Equation ([T]) implies that the interference between the 
c and cc running modes should give rise to sine and co- 
sine standing waves along the equator. The locations of 
the nodes and antinodes of and are automatically 
established by the random distribution of a large number 
of inhomogeneities in the silica sphere [1, ■ To visu- 
alize this effect, we have scanned a sharp fiber tip along 
the equator (i.e. in the (j) direction) of an already split 
fundamental WGM and have recorded spectra at each 
point (note that the radial coordinate of the tip is kept 
constant using a shear- force feedback [iB|). Figure [2^ 
shows that the intensities of the two peaks of a doublet 
undergo out of phase periodic modulations as a function 
of the tip location. A slight slope of the middle base line 
is attributed to a drift in the shear-force tip-sphere stabi- 
lization. At location (i), the tip is positioned in the node 
of the symmetric mode and the antinode of the antisym- 
metric mode. Thus, it induces loss in the E^ mode while 
it leaves E- nearly unaffected. Position (iii) shows the 
opposite counterpart of (i) whereas at position (ii) both 
modes are affected equally strongly. 

As shown in Fig. [^b, the three spectra reveal that 
in addition to a change in the intensity balance of the 
doublets, their splittings are also modified. Figure [Hi 
displays the modulation of the splitting about its initial 
value of 24 MHz shown by the dotted line. Interestingly, 
we find that the tip can not only increase the mode split- 
ting, but it can also decrease it. This is due to a destruc- 
tive interference between the field scattered by the tip 
and the field scattered by the inhomogeneities in the mi- 
crosphere that gave rise to the initial splitting. Figures[2^ 
and f provide further data on the increase and decrease 
of the mode splitting as the tip was scanned in the 9 di- 
rection for two different (p positions spaced by half of the 
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FIG. 2: a) The intensities of the two peaks in a doublet as 
a function of the tip position along the equator. The dotted 
line marks the slight intensity drift in the detection, b) Spec- 
tra recorded at positions (i-iii) in Fig. (a). The spectra are 
displaced vertically for clarity, c) Simultaneously measured 
spectra on PD2 and PD3. d) The recorded splitting corre- 
sponding to the data in Fig. (a). The solid curve shows a 
sinusoidal fit. e) and f) The variations of the mode splitting 
as a function of the tip in the 6 direction for positions (i) and 
(iii), respectively. The solid curves display fits according to 
the spatial mode function of the fundamental WGM. 



interference period along the equator. Finally, Fig. [2j; 
plots the resonance spectra recorded simultaneously on 
PD2 and PD3, i.e. via global scattering from the sphere 
and via the fiber tip. The different lineshapes on the two 
channels might seem unexpected at first. However, this 
effect shows that if the tip is placed in an antinode of E- 
or it efficiently extracts photons out of that mode, 
leading to a larger signal in the fiber tip and thus a lower 
intensity in the cavity mode. On the contrary, the mode 
that is less perturbed is stronger in the resonator and 
is nearly uncoupled to the fiber tip. It is evident that 
the mode that is coupled to the tip has experienced an 
additional broadening. 

We now show that the radiation properties of a sub- 
wavelength object such as its scattering rate are modified 
much in the same manner as those of the spontaneous 
emission of an atom. Our guiding thought is that many 
central features of CQED, including the modification of 
the mode density in a resonator, can be traced to the spa- 
tial character of the modes and should be thus shared by 
classical cavity electrodynamics. We first present a sim- 
ple treatment of the free-space Rayleigh scattering using 
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a semi-quantum electrodynamic (semi-QED) approach, 
where the material scatterer is treated classically while 
the field is quantized. Then we will discuss the modifica- 
tion of the scattering rate when the scatterer is coupled 
to a resonator. 

Let us assume that a freely propagating photon is 
incident on a subwavelength spherical scatterer of ra- 
dius a and refractive index n. We take the photon to 
be in a mode with volume Vk, frequency tOk and a 
linear polarization along the unit vector such that 
Efc = Skek{a\ -f Cfe) where a\. and are the usual cre- 
ation and annihilation operators and £k = yJfiUkj2eQVk- 
In the limit where the scatterer is considerably smaller 
than A, it can be described by a dipolar polarizability 
a [20| so that the induced dipole moment operator reads 
Pfc = eQa£k{a\ +ak)ek- Thus the interaction energy be- 
tween this dipole moment and another mode Ej becomes, 



(2) 



if we neglect the terms that do not conserve photon 
numbers. Here we have taken 



9k] 



/ '^^fVkVvac and have set Vj — Vuac for all vacuum 
modes j into which the incident beam is scattered. The 
system Hamiltonian becomes [2l| 



H = hujkaiak 



j 



E 



hgkjialcLj+a^ak), (3) 



leading to the Heisenberg equation of motion 



iak = LOkak + ^ 9kjaj. 



(4) 



The last term in Eq. (|4]) signifies the scattering of the 
incident field into all vacuum modes. Following the 
Weisskopf-Wigner formalism [22|, we find the rate 



(5) 



for this scattering event [2l|. Here we have restricted 
ourselves to tok = tUj for elastic scattering and have used 
the notation gn — —aLUk/2^/VkVvac- Now we can cal- 
culate the Rayleigh scattering cross section Q by 
considering the total power radiated by the scatterer ac- 
cording to /incCfi = Ti^k^R- Given that Lmc = 'huikc/Vk 
and a = Ana^ I "oTn I , one obtains the well-known relation 



I n2+2 I' 
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(6) 



We note that a rigorous quantum optical treatment of 
scattering is not frequently discussed in the literature [2^ 
and goes beyond the scope of this paper. However, the 
fact that (Jn can be derived via the Weisskopf-Wigner 
formalism using quantized fields provides a robust sup- 
port for the intuitive expectation that a modification of 



the mode density, for example in a resonator or in front 
of a mirror, could also lead to a change in the Rayleigh 
scattering rate. The corresponding Purcell effect offers a 
physical explanation for the question posed earlier. The 
rate with which energy is transferred from Ec to Ecc is 
enhanced by F and is given by 7]F where 77 is the ge- 
ometric factor determining the fraction of the solid an- 
gle subtended by the mode. Therefore, a fundamental 
WGM with Vm — 130 nrrv" (sphere diameter 30 nm) and 
Q — 10^ yields F ~ 10'*, compensating for the very small 
geometric acceptance of the order of 10"''. The influ- 
ence of the Purcell factor F also explains why reducing 
the cavity Q results in the disappearance of light in the 
counterpropagating mode as reported previously [l3| . 
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FIG. 3: a) The mode splitting measured as a function of 
sphere-tip separation in the r direction, b) Plot of the mea- 
sured splitting versus the additional tip-induced broadening. 
The solid curve displays the fit to a quadratic function. 

Having shown that the density of states plays a cen- 
tral role in the description of Rayleigh scattering, next 
we consider the coupling of two counterpropagating cav- 
ity modes Ec and Ecc via a single Rayleigh scatterer. 
The details of the calculations are pro vided in the sup- 
plementary materials of this paper [2ll |. Returning to a 
classical notation, we find 



(7) 



- ujc shows the 
is the mode 



E^ = (-i A - 70) E^ + Ko 
E+ = {-iA + 2ig - 70 - F) £;+ + kq. 

Here we have defined E = e*"*i?, A = 
detuning of the laser frequency lu, k — kqc 
excitation rate, and 270 denotes the unperturbed cavity 
linewidth. When dealing with Rayleigh scattering out of 
a cavity, we have to take into account the spatial variation 
/(r) of Ec and Ecc in the resonator mode. Going back 
to the definition of gkj and noting that Vj = Vk is the 
WGM volume Vm of the two modes, we thus obtain [2l| 



2g = -ap{r)ujc/Vm 
F = a^p{r)u;t/67Tc^Vm 



(8) 



for the mode splitting and broadening, respectively. 

A close scrutiny of Eq. ([7]) shows that if \2g\ is suffi- 
ciently large to overcome 70 and if 2| ^■2^2 1 ^ (2^0)'^ *° 



assure that \2g\ > F, a mode splitting is resolved. This is 
similar to the case of the strong coupling in CQED where 
the coherent exchange of energy between the cavity mode 
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and an atom leads to a mode splitting if the coupling co- 
efficient \g\ becomes larger than the cavity linewidth. As 
in CQED, the mode splitting \2g\ grows with decreasing 
Vm and increasing w, but the atomic dipole strength is 
replaced here by a, and the splitting is asymmetric. The 
antisymmetric mode remains unperturbed (see Eq. ([7])) 
whereas the symmetric mode undergoes an additional 
broadening given by T. This is a consequence of the fact 
that the phases of the new eigenmodes are self-adjusted 
so that the local scatterer is placed in a node (antin- 
ode) of the antisymmetric (symmetric) mode. This un- 
balanced splitting has been recently also observed for the 
coupling of a far-detuned cold atom ensemble to a high-Q 
ring cavity [3] ■ Indeed, in that case several million atoms 
detuned from their transition resonances also behave es- 
sentially as a large dielectric object. 

Equations ([8]) predict a linear dependence between the 
tip- induced splitting and line broadening if /(r) is varied. 
However, due to the finite spatial extent of the tip and the 
WGMs, a effectively grows as the tip enters the mode. 
This position dependence leads then to a quadratic re- 
lationship between the tip-induced broadening and split- 
ting. To investigate this effect, we have moved the tip 
in the radial direction and have recorded spectra at each 
location. Figure [3^ shows that as expected, the mode 
splitting becomes larger when the tip enters the evanes- 
cent field of the microsphere. Figure plots the tip- 
induced splitting versus the increase in the linewidth of 
confirming their quadratic dependence. 



To realize an ideal scenario for studying the interaction 
of a single well-defined scatterer with the fundamental 
mode of a microsphere, we have searched for spheres in 
which no mode splitting was observable in the beginning 
and it was created only when the tip was introduced. 
Figure |4] presents four spectra of a resonance as the tip 
was scanned in the 0-direction from outside the mode 
(a) to the maximum of the mode at the equator (d) at 
constant separation from the sphere surface. The 
mode is shifted in frequency by 13 MHz and broadened 
by about 6 MHz. Given that |25|/r = SX^/in^a ac- 
cording to Eqs. ([H]), the observed ratio of the splitting to 
broadening implies a radius of a ~ 140 nm for a spheri- 
cal Rayleigh scatterer. This is in very good quantitative 
agreement with the experimental parameters, consider- 
ing a typical value of 50-100 nm for the radii of curva- 
ture of SNOM tips and accounting for the overlap of the 
conical tip taper with the evanescent part of the mode. 

In conclusion, we have considered the phenomenon of 
Rayleigh scattering both in free space and in the presence 
of a resonator. By using a semi-QED approach, we have 
pointed out the roles of the modification of the density 
of states and of the Purcell effect in classical scattering. 
Our results demonstrate that although the introduction 
of a scatterer into a high-finesse resonator might be com- 
monly thought to introduce losses, it can mediate a co- 
herent coupling of the resonator modes and cause their 
consequent normal mode splitting. 
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FIG. 4: a-d) WGM spectra recorded at different 6 values of 
the tip location. At position a) the tip is nearly outside the 
spatial profile of the mode whereas at position d) it is in the 
mode maximum, i.e. the sphere equator. 



We are grateful for financial support from the Deutsche 
Forschungsgemeinschaft SP1113, the Swiss National 
Foundation (ETH group), NaFoG, Bcrhn (A. M.) and 
the Alexander von Humboldt Stiftung (L. de S. M.). 



[1] P. R. Berman, Cavity Quantum Electrodynamics (Aca- 
demic Press, 1994). 
[2] P. Domokos, M. GangI, and H. Ritsch, Opt. Commun. 

185, 115 (2000). 
[3] J. Klinner, M. Lindholdt, B. Nagorny, and A. Hem- 

merich, Phys. Rev. Lett. 96, 023002 (2006). 
[4] T. Aoki, et al., Nature 443, 671 (2006). 
[5] K. Hennessy, et al.. Nature 445, 896 (2007). 
[6] V. B. Braginsky, M. L. Gorodetsky, and V. S. Ilchenko, 

Phys. Lett. A 137, 393 (1989). 
[7] L. CoUot, V. Lefevre-Seguin, M. Brune, J. Raimond, and 

S. Haroche, Eur. Phys. Lett. 23, 327 (1993). 
[8] T. J. Kippenberg, S. M. Spillane, and K. J. Vahala, Opt. 

Lett. 27, 1669 (2002). 
[9] M. L. Gorodetsky, A. D. Pryamikov, and V. S. Ilchenko, 

J. Opt. Soc. Am. B 17, 1051 (2000). 



[10] D. S. Weiss, et al., Opt. Lett. 20, 1835 (1995). 

[11] J. Krupka, P. Blondy, D. Cros, P. Guillon, and R. G. 

Geyer, IEEE Trans. Micro. Theo. Tech. 44, 1097 (1996). 
[12] M. BorseUi, T. J. Johnson, and O. Painter, Opt. Exp. 

13, 1515 (2005). 
[13] K. Srinivasan and O. Painter, Phys. Rev. A 55, 023814 

(2007). 

[14] R. J. C. Spreeuw and J. P. Woerdman, Physica B 75, 96 
(1991). 

[15] A. Venugopalan, D. Kumar, and R. Ghosh, Pramana J. 

Phys. 40, 107 (1993). 
[16] P.-Y. Bourgeois and V. Giordano, IEEE Trans. Microw. 

Theo. Tech. 53, 3185 (2005). 
[17] S. Gotzinger, et al., Nano Lett. 6, 1151 (2006). 
[18] A. Mazzei, S. Gotzinger, L. de S. Menezes, V. Sandogh- 

dar, and O. Benson, Opt. Comm. 250, 428 (2005). 



5 



[19] S. Gotzinger, O. Benson, and V. Sandoghdar, Opt. Lett. 
27, 80 (2002). 

[20] D. Jackson, Classical Electrodynamics (Wiley and Sons, 
1999). 

[21] The EPAPS Supplementary material presents the details 
of the calculations leading to equations ((5} and (O. 



[22] P. W. Milonni, The Quantum Vacuum (Academic Press, 

New York, 1994). 
[23] B. J. Dalton, S. M. Barnette, and P. L. Knight, J. Mod. 

Opt. 46, 1107 (1999). 



